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I . INTRODUCTION 


The  feasibility  of  developing  an  optical  miss  distance  sensor  has 
been  investigated  recently.  Field  tests  and  a wind  tunnel  test  yielded 
inconclusive  results,  indicating  a need  for  a better  understanding  of 
the  optical  problems  involved.  The  work  reported  here  addresses  this 
need. 


The  miss  distance  sensor  can  be  regarded  as  a replacement  for  a 
tracer  round.  The  general  idea  is  that  an  optical  signal  is  sent  to 
the  projectile,  then  reflected  from  the  base  to  a detector.  The  light 
source  and  detector  are  both  located  near  the  gun;  consequently,  the 
light  will  travel  primarily  through  the  length  of  the  wake  of  the 
projectile.  For  the  analysis  of  the  system  it  is  important  to  know  how 
the  wake  affects  the  light  rays,  particularly  in  regard  to  the  strength 
of  the  signal  received  by  the  detector.  A typical  mean  density  distri- 
bution in  the  wake,  in  the  radial  direction,  shows  the  lowest  density 
on  the  axis.  From  this  it  can  be  deduced  that  the  wake  should  act 
like  a negative  lens,  tending  to  diverge  the  light  rays.  This  qualita- 
tive result,  first  described  in  Reference  1,  was  given  quantitative 
meaning  by  the  experiment  reported  there.  Although  the  method  of 
detecting  the  effect  is  not  the  same  in  this  experiment  as  in  the  miss 
distance  sensor,  the  experiment  is  important  because  it  provides 
accurate  results  in  a controlled  and  convenient  manner.  The  calcula- 
tions described  here  are  applicable  to  the  miss  distance  sensor  and  the 
experiment,  but  only  application  to  the  latter  is  discussed.  The 
former  will  be  considered  in  a separate  report. 

The  experiment  produces  end-view  shadowgraphs  of  the  projectile 
and  wake.  The  shadow  of  the  projectile  appears  magnified  because  it  is 
formed  by  the  diverging  light  from  a near  point  source  and  because  the 
wake  has  a negative  lens  effect.  The  experimental  set-up  is  sketched 
in  Figure  1.  The  light  source  is  pulsed  when  the  projectile  is  at  a 
distance* *  i from  the  photographic  film;  this  distance  is  determined 
accurately  by  means  of  a simultaneous  side  view  shadowgraph.  The  com- 
plete projectile  shadow  does  not  appear  in  the  end-view  shadowgraph 
because  of  the  hole  punctured  in  the  film.  A sample  result  (Case  # 1)** 
is  shown  in  Figure  2;  about  60  percent  of  the  circumference  is  discern- 
ible. 


1.  G.  D.  Kahl,  D . B.  Sleator,  and  D.  D . Shear , "Experiments  on  Wake 
Optical  Properties , 11  Ballistic  Research  Laboratories  Memorandum 
Report  No.  2496,  July  1975.  AD  B005619. 

* Definitions  of  quantities  are  given  in  the  List  of  Symbols  section. 

**  This  is  the  sixth  case  listed  on  p.  12  of  Reference  1;  also,  it  is 
Case  f listed  in  Table  I on  p.  13  of  this  report. 
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For  miss  distance  sensors  and  beam  rider  projectiles  or  missiles 
the  wake  through  which  the  radiation  travels  is  turbulent,  except 
possibly  for  a short  length  near  the  base.  Considerable  work  has  been 
done  on  propagation  of  electromagnetic  radiation  through  a turbulent 
medium  because  of  applications  to  the  atmosphere;  however,  the  approxi- 
mations appropriate  to  that  work  will  generally  not  apply  to  a turbulent 
wake.  An  important  question  is  whether  geometrical  or  physical  optics 
should  be  used.  The  answer  depends  on  the  optical  system  (whether  it  is 
diffraction  limited  or  not),  the  size  and  shape  of  the  important  turbu- 
lent eddies  (turbulence  scale  lengths),  and  the  density  fluctuations. 

In  this  report,  rather  than  justify  the  geometrical  optics  approximation 
a priori,  we  demonstrate  that  its  predictions  compare  favorably  with 
experimental  results  of  Reference  1.  The  additional  question  arises  as 
to  the  importance  of  scattering  which  can  also  be  treated  within  the 
geometrical  optics  approximation.  We  calculate  only  the  refractive 
effects  due  to  mean  density  variations  in  the  wake;  the  results  of 
Reference  1 are  used  to  justify  neglecting  scattering. 

The  calculation  requires  that  the  index  of  refraction  be  specified; 
therefore,  the  density,  using  the  Dale-Gladstone  law.  For  a turbulent 
wake  the  density,  p,  is  expressed  as 

P = P + p' 

where  p is  the  mean  density  and  p'  is  the  fluctuating  part  of  p.  The 
optical  effects  of  p*  would  be  to  scatter  or  diffuse  the  light;  varia- 
tions in  p will  produce  refraction,  the  only  effect  considered  here. 

An  axisymmetric  wake  is  assumed  in  this  work.  Actually  the  wake 
is  never  axisymmetric,  even  for  a projectile  at  zero  angle  of  attack. 

The  near  wake  is  usually  close  to  being  axisymmetric;  but  the  far  wake, 
i.e.,  many  calibers  from  the  body,  can  develop  oscillations  with 
amplitudes  of  the  order  of  the  wake  diameter2.  The  complications  of  a 
statistical  treatment,  even  for  p,  are  not  warranted  for  this  study, 
and  the  comparisons  with  results  from  Reference  1 justify  this  approach. 


II.  EQUATIONS 

Figure  3 shows  a diagram  of  the  coordinate  system,  the  origin 
being  in  the  center  of  the  projectile  base.  The  ray  tracing  equations 
are  derived  from  Fermat fs  principle,  which  states  that  a ray  of  light 
will  travel  between  two  points  along  the  path  for  which  the  time  of 
travel  is  an  extremum.  (All  lengths  have  been  non-dimensionalized  by 


2.  C.  H.  Murphy  and  E.  R . Dickinson 3 " Growth  of  the  Turbulent  Wake 
Behind  a Supersonic  Sphere 3 n AIAA  Journal 3 Vol.  13  pp.  339-342 ^ 
February  1963 . Also  see  BEL  Memorandum  Report  No.  1388 . AD  280216. 
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D,  the  diameter  of  the  projectile  base;  dimensional  lengths  will  be 
designated  with  a bar.)  For  this  variational  problem  the  resulting 
Euler  equations3  are 


3F/3x  - d(3F/3x)/dz  = 0 
3F/3y  - d(3F/3y)/dz  = 0 

where 


(1) 


x = dx/dz,  y = dy/dz 

F (z , x,  y,  x,  y)  = n(x,  y,  z)  (x2  + y2  + l)'2  (2) 

For  an  axisymmetric  index  of  refraction  field,  n(x,  y,  z)  = 
n(r,  z) , the  equations  can  be  simplified  if  cylindrical  coordinates 
are  used: 

x = r cos  0,  y = r sin  0,  z = z 

In  this  formulation  0 is  an  ignorable  coordinate,  in  the  language 
of  dynamics.  Although  n is  axisymmetric  the  light  rays  will  be  three- 
dimensional  curves  in  general;  these  are  called  skew  rays.  A ray 
remains  in  a plane  only  if  the  initial  direction  is  in  a plane  through 
the  axis;  these  are  called  planar  rays.  It  can  be  shown4  that  a skew 
ray  can  be  treated  as  a planar  ray  if  the  index  n is  replaced  by 

m = [n2  - (h2/r2)  ~\  2 , 

where  h is  a constant.  Specifically,  the  Euler  equation  for  r(z)  is 


d2r/dz2  = [1  + (dr/dz)2]  [3m/3r  - (dr/dz)  (3m/3z)]/m  (3) 

with  0 (z)  determined  from  the  quadrature 

d0/dz  = (h/r)  [1  + (dr/dz) 2]  2/ (n2r2  - h2) 2 . (4) 

The  quantity  h,  constant  along  a light  ray,  can  be  related  to  the 
location  and  direction  of  the  ray  at  the  initial  point,  z = zq;  thus 

h = noro2(d0/dz)Q/[l  + (dr/dz)Q2  + rQ2  (d0/dz)Q2]'2  . 


3.  M.  Born  and  E.  Wolf , Princvp les^  of  Optics , Pergamon  Press , New 
York , 1959 , p.  717. 

4.  R.  K.  Luneburg , Mathematical  Theori)  of  Optics,  University  of 
California  Press , Berkeley  and  Los  Angeles , 1964,  p.  191. 
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The  Euler  equations  simplify  further  for  the  set  of  planar  rays;  for 
these  h = 0 and  d0/dz  = 0 from  Eq.  (4). 

The  index  of  refraction,  n(z,  r)  , together  with  its  partial 
derivatives,  is  provided  as  input  to  the  calculation. 


III.  DENSITY  DISTRIBUTION  MODEL 

The  index  of  refraction  in  the  flow  field  is  related  linearly  to 
the  density  of  the  gas: 

n = noo  + (n«,  - [Cp/pJ  - 1]  • (5) 

In  this  form  of  the  Dale-Gladstone  law,  the  wave  length  of  light  enters 
through  n , the  ambient  index  of  refraction.  A crucial  part  of  the 

problem  is  to  obtain  accurate  data  for  the  wake  density.  A semi- 
empirical  model  developed  by  Finson^  for  the  wake  density  distribution 
is  employed  in  this  work.  It  has  proven  satisfactory  for  the  cases 
studied  thus  far. 

The  key  to  Finson's  wake  model  is  the  specification  of  the  velocity 
on  the  axis.  Beyond  the  rear  stagnation  point,  the  mean  velocity  on  the 
axis  can  be  specified  as  a nearly  universal  function  of  downstream 
distance,  z.  This  "universal”  function  is  shown  in  Figure  2 of  Refer- 
ence 5,  based  on  the  results  of  nine  experimental  investigations. 
Additional  operations  are  prescribed  for  the  near  wake.  Conservation 
of  energy  can  then  be  used  to  relate  the  velocity  and  temperature  on 
the  axis.  Approximate  radial  profiles  can  be  specified  for  the  temper- 
ature from  similarity  arguments.  The  density  field  is  then  obtained 
with  the  perfect  gas  law,  with  the  additional  approximation  that  the 
pressure  is  constant,  equal  to  the  ambient  value.  The  formulas  used  to 
determine  density  and  density  gradients  are  presented  in  Appendix  A. 

The  Finson  model  is  expected  to  be  applicable  for  < 5,  for  relatively 

high  Reynolds  numbers  (104  - 105 6) , and  for  a wide  variety  of  body 
shapes . 


IV.  NUMERICAL  SOLUTION 

Refraction  ahead  of  the  projectile  base  is  neglected,  i.e.,  the 
light  rays  travel  undeviated  from  z = 0 to  the  film  with  the  direction 
attained  at  z = 0.  This  assumption  is  justified  because  (1)  the 


5.  M.  L.  Finson 3 "Density  Distributions  in  Supersonic  Projectile 

Wakes3"  PSI  TR-663  Physical  Sciences 3 Inc.3  Woburn 3 Massachusetts 3 
Ballistic  Research  Laboratories  Contract  Purchase  Order  No. 

DAAD05- 76-M-D1 11 3 October  1976. 
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distance  traversed  here  by  the  light  ray  is  small  compared  to  the  length 
of  the  wake  traveled,  and  (2)  the  expansion  at  the  edge  of  the  projec- 
tile base  and  the  bow  shock  wave  produce  opposing  deviations. 

The  objective  of  the  ray  tracing  calculation  is  to  determine  the 
shadow  of  the  projectile  on  the  film  ahead  of  it.  The  rays  are  deter- 
mined by  solving  Eqs.  (3)  and  (4),  with  initial  conditions  specified 
by  the  position  of  the  light  source  and  the  slope  of  a ray  emanating 
from  it.  A point  source  on  the  x-axis  is  assumed  with  coordinates  # 

z = z,r  = r,0  = O.  The  slope  of  a ray  is  identified  by  r and  0 
o O 0 0 

(dr/dz  and  d0/dz  at  z = z ) . The  numerical  integration,  for  given  fQ 
and  0 , is  terminated  at  z = 0.  The  ray,  for  which  r = rb,  0 = 0b  at 
z = 0,  is  extended  as  a straight  line  in  the  direction  r^,  0^;  and  its 
intersection  with  the  screen  occurs  at 

xs  = ^ rb  ®b  Sin  ®b  + (rb  " £ rb}  C°S  ®b 

(6) 

= (rb  ' £ V Sin  °b  " 1 rb  °b  C0S  9b  ’ 

where  l is  distance  from  the  base  of  the  projectile  to  the  screen. 

The  computations  were  carried  out  on  a UNIVAC  1108  and  a CDC 
Cyber  76.  The  integration  interval  was  maintained  sufficiently  small 
in  the  integration  of  Eqs.  (3)  and  (4)  so  that  the  solution  was  correct 
at  least  to  four  figures. 

All  but  one  of  the  cases  treated  here  apply  to  a caliber  30  bullet. 
Although  the  base  diameter  is  not  quite  equal  to  the  full  cross- 
sectional  diameter  of  the  bullet  because  of  a slight  boattail,  we  take 
it  equal  to  the  latter.  We  calculate  the  "shadow"  of  the  bullet  as  the 
trace  on  the  film  of  those  rays  which  touch  the  base  at  its  circum- 
ference; i.e.,  rfe  = 0.5.  We  call  these  the  "grazing  rays"  and  denote 

them  by  r ,0  . A systematic  procedure  for  determining  grazing  rays 

3 og’  og 

is  described  in  Appendix  B. 

Our  choice  of  shadow-bounding  rays  is  an  approximation;  some  rays 
touching  the  base  at  r^  = 0.5  are  blocked  by  the  bullet  and  do  not  truly 

constitute  grazing  rays.  This  effect  is  discussed  in  Section  V;  it  is 
concluded  that  the  shadow  boundary  is  not  seriously  affected  by  this 
approximation . 
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V.  PLANAR  RAYS 


Planar  rays  (mentioned  in  Section  II)  are  particularly  useful; 
they  simplify  ray  tracing,  but  results  are  still  informative.  Here 
0=0=0;  all  light  rays  are  in  the  x,  z plane  and  are  characterized 
by  a single  parameter,  r , for  a fixed  light  source  and  screen.  The 

grazing  rays,  r for  x^  = ± 0.5,  are  easily  found  from  the  computed  x^ 

vs  r relationship.  These  rays  produce  the  two  boundary  points  of  the 

projectile  shadow  on  the  x-axis  if  there  is  no  caustic.  Since  the 

shadow  of  the  projectile  is  approximately  a circle  on  every  shadowgraph, 

the  two  bounding  points  of  light  penetration  provide  an  indication  of 

the  size  and  displacement  of  the  shadow.  For  these  rays  the  computed 

magnification  of  the  projectile  cross-section  is  defined  as 

Ax  ==  x - x ,,,  x and  x - being  the  coordinates  on  the  screen  of 
mag  sn  sf  sn  sf 

the  fnear!  and  !farf  planar  grazing  rays,  respectively.  The  'near* 
grazing  ray  is  the  one  on  the  same  side  of  the  shadow  as  the  light 
source. 

For  Case  # 1 (Figure  2)  the  two  planar  grazing  rays  correspond  to 
(1)  r = 0.00941  (near  ray)  and  (2)  r^  = 0.01101  (far  ray).  Figure  4 

shows  these  rays;  the  z scale  was  compressed  in  order  to  exhibit  the 
full  ray  conveniently.  The  far  grazing  ray  changes  curvature  near  the 
z-axis  and  is  concave  downward  when  it  strikes  the  projectile  base. 

In  Table  I*  we  present  a comparison  of  measured  and  computed  magni- 
fications for  nine  cases.  The  first  eight  are  those  listed  on  p.  12  of 
Reference  1.  Rounds  a through  h are  caliber  30  bullets;  parameters  are 
D = 0.782  cm,  r^  = 5.754  projectile  diameters,  n^  = 1.0002780,  = 

0.316,  M = 2.41,  y = 1.4,  and  T /T  = 1.10.  Round  i is  a spherical 

3 CO  7 1 9 W 00 

projectile;  the  "base"  is  taken  to  be  the  area  of  intersection  with  a 

plane  (normal  to_  the  trajectory)  through  the  center  of  the  sphere. 

Parameters  are  D = 0.800  cm,  r = 2.063,  n = 1.0002780,  C = 0.968, 

o 00  D 

M = 2.81,  y = 1.4,  and  T /T  = 1.00. 

CO  YV  00 

With  the  exception  of  Case  h,  which  appears  to  be  out  of  line  in 
the  measurement  column,  the  experimental  and  theoretical  values  agree 
to  within  4*2%.  The  close  agreement  is  a little  surprising  considering 
the  approximation  involved  in  the  model,  especially  for  the  wake 
density.  The  consistency  of  the  data  and  calculations  is  taken  as 
justification  of  the  model. 


* The  dashed  lines  separate  the  three  values  of  z + L 
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Table  I. 

Measured 

and  Computed  Magnifications 

Case 

z + £ 

£ 

Meas.  Magnif. 

Ax 

mag 

0 

a 

259.9 

13.54 

1.140 

1.108 

b 

259.9 

15.17 

1.175 

1.121 

c 

646.4 

14.68 

1.150 

1.119 

d 

646.4 

26.34 

1.230 

1.213 

e 

646 . 4 

30.50 

1.260 

1.246 

£ 

646.4 

52.85 

1.410 

1.420 

g 

646.4 

75.13 

1.640 

1.603 

h 

646 . 4 

90.95 

1.630 

1.721 

i 

m-z 

31.84 

1.366 

1.373 

A parameter  study  was  carried  out  corresponding  to  an  experiment 
where  the  light  source  and  screen  are  fixed,  but  the  timing  of  the 
light  pulse  is  varied;  thus  rQ  and  z^  + £ are  fixed.  The  values  are 

those  of  Cases  c through  h in  Table  I.  Figure  5 shows  how  the  magni- 
fication, Axmag>  changes  with  location  of  the  projectile  and  illustrates 

strikingly  the  lens  effect  of  the  wake  by  comparing  the  wake  flow  and 
no  flow  magnifications.  It  is  noteworthy  that  the  wake  flow  ^xmag 

varies  linearly  with  £ over  a wide  range  of  £. 

Part  of  the  light  striking  the  base  at  its  edge  on  the  near  side 
is  blocked  by  the  projectile  so  that  a portion  of  the  edge  of  the 
shadow  is  formed  by  rays  other  than  those  passing  through  r^  = 0.5. 

An  estimate  based  on  planar  rays  can  be  made  of  the  deviation  of  our 

computed  shadow  from  that  formed  by  the  true  grazing  rays  on  the  near 

side.  We  assume  that  the  projectile  is  a cylinder  of  length  k having 

a blunt  nose  so  that  the  actual  near  grazing  ray  (planar)  passes  through 

r = 0.5  at  the  front  end  of  the  cylinder.  Then  the  corrected  shadow 

coordinate  x = (x  ).  can  be  obtained  in  two  steps: 

s v sn^i  r 

(1)  find  r . , the  value  of  r for  which 
v J oi’  o 
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0.5 


and  designate  the  corresponding  r^  by 

(2)  evaluate  (x  ).  from 
sn  1 

(x  ).  = r . - (£/k)  (r  . - 0.5) 
v siri  bi  J bi 

For  a length  k = 5,  the  "corrected”  magnifications  exceed  the 
original  values  by  no  more  than  1%%  in  Cases  d through  h in  Table  I. 
This  result  suggests  that  neglect  of  blocked  rays  does  not  produce 
serious  errors. 

An  inspection  of  the  data  was  made  to  determine  whether  any  non- 
blocked  light  rays  fell  inside  the  trace  of  the  grazing  rays  on  the 
screen.  This  occurrence  would  denote  the  presence  of  a caustic.  No 
caustic  was  found  in  any  of  the  cases  studied. 


VI.  SHADOW  OF  PROJECTILE 

The  shadow  of  the  projectile  is  bounded  by  the  locus  of  the 
incident  grazing  rays  on  the  screen.  Figure  6 shows  the  calculated 
shadow  for  Case  No.  1.  Also  shown  are  the  location  of  the  projectile 
in  the  x,  y plane  and  the  no- flow  image  produced  by  the  point  source. 
The  latter  is  a circle  given  by 

[xs  + (ro£/zo)]2  + ys2  = (1/4)  [1  + (Vzq)]2  , (7) 

where  the  subscript  s denotes  the  coordinate  value  on  the  screen.  The 
diagram  shows  that  the  wake  has  appreciably  enlarged  the  shadow  beyond 
the  simple  geometrical  projection,  or  no-flow  image.  The  image  is 
symmetrical  about  the  x-axis  because  the  point  source  is  located  on  the 
axis.  Although  it  is  not  a circle,  it  is  only  slightly  eccentric,  the 
height  being  larger  than  the  width  by  approximately  3^%. 

Figure  6 provides  the  theoretical  picture  to  be  compared  with  the 
shadowgraph  of  Figure  2.  The  projectile  shadow  in  Figure  2 does  not 
appear  clearly  symmetrical  about  any  axis.  This  lack  of  symmetry  is 
probably  caused  by  an  unsymmetrical  wake  and  the  yawing  of  the  bullet. 
Other  shadowgraphs  do,  however,  exhibit  more  pronounced  symmetry. 

The  calculation  predicts  two  effects  of  the  wake;  namely, 

(1)  magnification  of  the  object,  discussed  in  Section  V,  and  (2)  dis- 
placement of  its  image  to  the  left.  The  displacement  is  more  difficult 
to  measure  because  it  additionally  requires  locating  the  line  of  fire. 
Displacement  was  measured  for  only  one  of  the  rounds  listed  in  Table  I 
of  Section  V;  namely,  for  the  spherical  projectile,  Round  i.  The 
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measured  displacement  of  the  center  of  the  circle  was  0.509  calibers 
while  the  computed  value  was  0.525  diameters. 

An  additional  comparison  of  measured  and  computed  horizontal 
magnifications  was  made;  in  this  case  the  projectile  was  a sphere. 

The  grazing  rays  were  taken  to  be  those  which  intersected  the  circle 
formed  by  the  sphere  and  the  plane  (normal  to  the  trajectory)  through 
its  center.  The  two  results  disagreed  by  less  than  one  percent. 


VII.  CONCLUDING  REMARKS 

On  the  basis  of  results  reported  here  and  in  Reference  1,  it 
appears  feasible  to  make  optical  measurements  in  ballistic  experiments 
in  which  light  traverses  the  wake  of  a projectile.  Furthermore,  it 
seems  possible  to  predict  the  paths  of  light  rays  with  reasonable 
accuracy  on  the  basis  of  geometrical  optics  and  theoretical  models  for 
wake  flow.  The  model  presented  here  is  verified  by  experimental 
results  and,  therefore,  can  be  used  with  confidence  in  performing 
design  calculations  in  the  development  of  miss  distance  sensors  and 
beam  rider  projectiles  and  missiles. 
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Figure  1.  Physical  Setup  of  Experiment 


Figure  2.  End-On  Shadowgraph  of  Projectile  --  Case  No.  1 
(D  = 0.782  cm,  H = 41  cm,  zq  = 464  cm, 

r = 4.5  cm,  M = 2.41,  C = 0.316) 
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Figure  3.  Coordinate  System 


D = 0.00782  m r0  1 5.754 


SOURCE 


Figure  4.  Planar  Grazing  Rays  --  Case  No.  1 
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Figure  5.  Magnification  of  Projectile  Base  as  a Function  of  Location 
of  Projectile  (Light  Source  and  Screen  Locations  Fixed)  -- 
Case  No.  1 
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Figure  6.  Calculated  End-On  Shadow  of  Projectile  --  Case  No.  1 
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LIST  OF  SYMBOLS 


Dimensional  lengths  are  designated  by  a super  bar;  otherwise 
lengths,  including  coordinates,  are  non-dimensionalized  by  D. 


D 

F 

h 

H 

k 

& 

m 

M 

< 

n 

n 


projectile  drag  coefficient 
diameter  of  projectile  base 

integrand  in  extremum  integral  --  Eqs.  (1)  and  (2) 
constant  on  a three-dimensional  light  ray  --  Eq.  (3)  et  seq. 
total  enthalpy  of  air 

length  of  cylindrical  portion  of  projectile 

distance  between  base  of  projectile  and  downrange  screen 

equivalent  refractive  index  in  planar  ray  equation  for  non- 
planar  rays  (Section  II) 

flight  Mach  number  of  projectile 

index  of  refraction  of  air  in  flow  field 
ambient  index  of  refraction  of  air  (=  1.0002780) 


r,  6 

r , 0 
og*  og 


w 


cylindrical  polar  coordinates  (0  in  radians)  --  Figure  3 
initial  values  of  r and  0 for  grazing  rays 

temperature 

an  average  temperature  of  projectile  wall 
ambient  temperature 


u 


^ , y 

s s 


X , X 

sn 


X,  y, 


sf 

z 


mean  axial  velocity  component  relative  to  projectile 
x,  y of  a point  on  the  film 

values  x^  for  near  and  far  planar  grazing  rays,  respectively 
rectangular  coordinates  --  Figure  3 
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LIST  OF  SYMBOLS  (Continued) 


^Xmag  (=  xsn  - xsf)  magnification  of  projectile  base 

Y ratio  of  specific  heats  of  air  (=  1.4) 

p flow  field  density  of  air 

Subscripts 

0 z = z --  plane  of  light  source  (initial  value  of  variables) 

b z = 0 — plane  of  projectile  base 

c£  r = 0 --  center  line 

g grazing  rays  (for  which  r^  = 0.5) 

1 planar  rays  grazing  the  front  end  of  cylinder 

s plane  of  film  (screen)  --  z = - £ 

00  free  stream,  or  ambient,  conditions 


Superscripts 

= d/dz 

— dimensional  length 
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APPENDIX  A:  REFRACTIVE  INDEX  IN  THE  WAKE 


The  following  quantities  are  assumed  known:  M^,  C^,  y,  n^,  and 

T /T  . The  index  of  refraction,  n(r,z),  in  the  flow  field  behind  the 

w 00  v 

projectile,  described  in  Reference  5,  is  obtained  from  the  formulas 
which  now  follow,  evaluated  in  sequence. 


The  formulae  for  computing  the  velocity  deficit  are  given  in 
Eqs.  (A. 1) - (A. 13) . The  velocity  deficit,  fi,  expressed  as  u = u^  + 

where  is  an  analytical  approximation  to  the  empirical  curve  for 
Au^/u^  in  Figure  2 of  Reference  5,  and  u^  is  an  approximation  to  the 

curve  in  Figure  3 of  Reference  5. 


u = (u  - u „)/u  --  definition,  where  u is  mean  axial  velocity 

00  c£  OO  . . 

component  relative  to  projectile 


z = z (1.1284/Cd2)  , 


dz/dz  = 0.88623  CD 2 


If  0 < z < 1 , 


u:  = 0.5  , 


du^dz  = 0 


(A.  1) 
(A.  2) 


If  1 < z < 13.0  , 


= % exp  [-0.072662  (In  z)2] 


du^/dz  = -0.14532  (In  z)/z 


(A.  3) 


If  13.0  < z < 1,000.0  , 


If  1,000.0  < z 


* ->  •zcti  *-0.79073 

= 2.3561  z 


du^/dz  = -0.79073  u^/z 


*-2/3 


= z 


* , * * . * 
du^/dz  = - (2/3)  u^/z 


(A. 4) 


(A. 5) 


q = z - 
If  q < 0.0  , 


U2  = 0.50  , 


, * . , 1c  _ _ 

du^/dz  = 0.0 


(A.  6) 
(A.  7) 


If  0.0  < q < 0.060  , 


u2  = 0.50  - (13/3)  q 


du2/dz  = - (13/3)  dz/dz 


(A.  8) 
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(A.  9) 


If  0.06  < q < 0.340 


If  0.340  < q < 0.671 


u2  = 0.34083  - 1.8761  q + 3.2618  q2 
du^/dz  = [-  1.8761  + 6.5236  q]  dz/dz 
u2  = 0.12109  - 0.12085  q 


du2/dz  = - 0.12085  dz/dz 


(A. 10) 


If  0.671  < q < 1.321 


u2  = 0.081292  - 0.061539  q 
du^/dz  = - 0.061539  dz/dz 


(A. 11) 


If  1.321  < q , 


u2  = 0 , 


du2/dz  = 0 


* * * 

U = U!  + U2  5 


du/dz  = du^/dz  + du2/dz 


(A. 12) 
(A. 13) 


The  relation  between  velocity  deficit  and  temperature  is  developed 
in  Eqs.  (A. 14)- (A. 18) . Eq.  (A. 18)  is  Eq.  (2)  of  Reference  5.  The 
centerline  deficit  in  total  enthalpy,  AH  in  Eq.  (A. 17),  is  small  for 
the  applications  of  this  model  to  projectile  flight.  We  have  derived 
a different,  presumably  more  accurate,  expression  for  it;  viz, 

if 

AH  = 0.0136  R 2 (2d/D)  [1  - (T  /T  ) N _1]  , 

W ^Vv<»  ooJ 

where  (T  is  an  effective  heat  transfer  body  length,  approximately  equal 
to,  but  less  than,  the  actual  body  length. 

* 

AH  = (H^  - H p/H^  --  definition  where  H is  total  enthalpy  and 


h is  enthalpy:  H = h + u2/2. 


T = T /T  --  definition 
c£  °o 


If  z < 0. 5 


R2  = 26.514  Cp-2/3  r2 


R = R (r  = %) 
w *•  2J 


R 2 = 6.6285  f 
w D 


-2/3 


(A. 14) 


dR  2/dz  = 0 
w 
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If  0.5  < z , R2  = 16.703  r2  (Cp  z)~2^ 


R 2 = 4.1758  (Cn  z)"2/3  (A. 15) 

w D 

dR  2/dz  = -(2/3)  R 2/z 
w ' J w 

N00=1+js(Y-1)  M^2  (A.  16) 

AH  = 4 x 10'3  Rw2  [1  - (T/TJ  N^'1]  (A.  17) 

T = 1 + % (Y  - 1)  Mro2  (2  - u)  u - NJM  (A. 18) 


The  a that  is  obtained  from  Eq.  (A. 19)  is  the  n of  Eq.  (5)  in 
Reference  5.  This  is  a density- transformed  radius,  depending  both  on 
r and  z.  Eq.  (A. 22)  results  from  the  assumption  of  a Gaussian  radial 
profile  of  temperature  with  respect  to  a,  the  perfect  gas  law,  and  the 
assumption  of  constant  pressure. 

Now  solve  the  following  functional  equation  for  a: 

F (a)  = a + 2.941  (T  - 1)  [1  - exp  (-0.34  a)]  - R2  = 0 (A. 19) 

This  equation  is  solved  by  the  method  of  false  position;  the  following 
two  initial  guesses  a and  a have  proven  successful  in  generating  a 

convergent  sequence  of  approximations: 


Oj  = R2  - 2.941  (T  - 1)  , 

o2  = R2/T 

(A. 20) 

E = exp  (-0.34  a) 

(A. 21) 

p/Poo  = 1/ [1  + (T  - 1)  E] 

(A. 22) 

n = noo + (A*  - n tcp/pj 

- 1] 

(A. 23) 

3 (p/pj/3r  = 0.68  (R2/r) 

(P/P J3  (T  - 1)  E 

(A. 24) 

d (AH) / dz  = 4 x 10-3  [1  - 

(T  /T  ) N _1]  dR  2/dz 

v yj*  go'  oo  j yj 

(A. 25) 

Q = (y  - 1)  (1  - u)  du/dz 

(A. 26) 

* 

* 

dT/dz  = 1.1284  (E,  2 Q - N 

D « 

_ d (AH) / dz 

(A. 27) 
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9 (p/pj/3z  = - (p/p J2  E [1  + (T  - 1)  (1  - E)  (p/pj]  dT/dz  (A. 28) 
9n/9r  = (n^  - 1)  9 (p/poo)/9r  (A. 29) 
9n/9z  = (na  - 1)  9 (p/pco)/9z  (A. 30) 

h = h (r,  z , r,  0)=nr20/[l  + r2  + r 20  2]2  , 

where  nQ  is  evaluated  from  Eq.  (A. 23)  for  rQ,  zq  the  very  first  time 
this  sequence  of  formulas  is  used.  (See  Eq.  (4),  et  seq.) 

m = [n2  - (h/r)2]  2 (A. 31) 
9m/ 9 z = (n/m)  9n/9z  (A. 32) 
9m/9r  = (1/m)  [n  9n/9r  + (h2/r3)]  (A. 33) 


26 


APPENDIX  B:  DETERMINATION  OF  GRAZING  RAYS 


For  a given  firing,  r^  = r^  (fQ,  0 ),  and  the  locus  of  grazing 

rays  in  the  r , 0 plane  must  be  found  numerically.  The  usual 

procedure  is  to  take  a fixed  value,  say,  of  0q,  then  compute  r^  for  a 

range  of  r and  pick  that  value  of  fQ  for  which  r^  = 0.5.  Reasonable 

estimates  for  r and  0 can  be  found  beforehand,  as  shown  below, 
og  og 


A very  useful  aid  is  a plot  of  the  r , grazing  ray  locus  for 

the  no-flow  condition,  namely  for  the  straight  lines  between  the  light 
source  and  the  boundary  of  the  projectile  base.  *The  equation  is 


02 

og 


[(1/4)  - (r 


r z)^]/(rz)2 
og  oJ  J K o oJ 


This  curve  is  shown  in  Figure  B.l  for  Case  No.  1. 

The  appropriate  r , 0q^  curve  can  be  obtained  easily,  e.g.,  for  a 

succession  of  points  moving  gradually  away  from  the  planar  values  of 

r . For  a given  0 the  range  of  r to  be  searched  can  be  estimated  in 
Og  O 0 

advance  (Figure  B.l).  As  an  example,  in  Figure  B.l,  for  0q  = 0.2  x 10 

the  two  dashed  horizontal  segments  represent  reasonable  search  ranges 

for  r . In  other  areas  it  will  be  more  appropriate  to  fix  r and 
og  og 

search  for  0 , as  for  example,  at  r = 0.0100. 

og*  r og 
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Figure  B.l.  Locus  of  Grazing  Rays  in  fQ,  0q  Plane  --  Case  No.  1 
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